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Abstract 

We present the results of the analytical evaluation of the massless four-loop 
0(a^) correction to the correlator of the quark scalar currents and the Higgs 
decay rate into hadrons. In numerical form we found (in MS scheme) that 
T(H -> 56) = ^M H m 2 b (M H ) [1 + 5.Q67a s /n + (35.94 - 1.359n/) (a s /ir) 2 

+ (164.139 — 25.771nj + 0.2h9rvj) (a s /7r) 3 where rif is the number of quark 
favours and a s = a s (M H ). 
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1 Introduction 



In the Standard Model (SM) one physical scalar Higgs boson is present as a remnant 
of the mechanism of mass generation. Particularly interesting for the observation of the 
Higgs boson with an intermediate mass Mjj < 2M^/ is the dominant decay channel into a 
bottom pair H —> bb. ( Standard Model properties of the Higgs boson have been discussed 
in many reviews; see for example |l|, £§). The partial width T(H — > bb) is significantly 
affected by QCD radiative corrections. First order a s corrections including the full m& 
dependence were studied by several groups || |, [|, ||, [/J. Second order corrections were 
calculated in the limit ml <C Mjj. Apart from the trivial overall factor ml due to the 
Yukawa coupling, corrections were obtained for otherwise massless quarks in Refs. || and 
for a nonvanishing mass of the virtual top quark in Ref. (both results are confirmed 
in Ref. iHj). Subleading corrections in the ml/Mjj expansion were found in Ref. [fLlfl 
and also confirmed in Ref. [fLOfl . In the latter work an additional quasi- massless (and 
numerically important) contributions of order a 2 have been identified and elaborated. 
They come from so-called singlet diagrams with a non-decoupling top quark loop inside 
(similar effects for the decay of the Z-boson to quarks were first discovered in pi). The 



results of Ref. |L0[ have been confirmed and extended by computing extra terms in the 



expansion in M^/m 2 in Ref. [[13 

In this work we compute the next-next-to-leading massless correction of order a z s to 
T(H —> hadrons). 



2 Preliminaries 

We start with considering the two-point correlators 

Tl s ff ,(Q 2 ) = {A-nfi J dxe* x {0\ T[ J S f (x) jf,(Q) ] |0>. (1) 

Here Q 2 = —q 2 , J| = ^ff^f is the scalar current for quarks with flavour / and mass mf, 
which are coupled to the scalar Higgs bosons. The total hadronic decay rate of a scalar 
Higgs boson is naturally expressed in terms of the absorptive part of the correlator ([!]) 

R s ffl (s)^^- s hnU s ffl (-s-ie), 



as follows §|,§: 

T(H -> hadrons) = -^^M H Y,m f m r R b ff ,(s = M 2 H ). (2) 
4V27T jj, 

Assuming that the Higgs boson mass Mh is less then 2m t and much larger than mj, we 
shall treat correlators (|1|) within the massless five-quark QCD. (In fact, even in the limit 
of heavy top the Higgs boson coupling with the current Jf still contributes to Th', we 



shall ignore these effects in our work. For a recent discussion see e.g. Ref. [14].) Under 
such an assumption the nondiagonal correlators n^j, with / ^ f vanish identically while 
the diagonal correlators get flavour independent, that is 

n^(Q 2 ) = n s (Q 2 ) and R s ff (s) = R s (s). 
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The renormalization mode of the polarization operator II s (Q 2 ) reads (see, e.g. Ref. 

n s (Q 2 ,a„/i 2 ) =ZfQ 2 + Z 2 m Yil{Q\a% (3) 

where a s = a s /ir = g 2 /(47r 2 ), g is the strong coupling constant; Z m is the quark mass 
renormalization constant. Within the MS scheme |16| 



1<j<» 



and the coefficients (Z^ s ) are just numbers, with D = 4 — 2e standing for the space-time 

dimension. As a result we arrive at the following renormalization group (RG) equation 
for the polarization operator II s (Q 2 ) 

[fW 2 + asf3{as) ^k + 2lm{as) ) nS = 7 " S(as) (5) 

2 

or, equivalently, (L Q = In ^ ) 

3 -n s = 7 f (a s ) - f 2 7m (a s ) + (3(a s )a s ^-) II s . (6) 



oLq 1 \ da S/ 

Here the anomalous dimensions 7g S (a s ) and 7m (a s ), and the /3-function (3(a s ) are defined 
in the usual way 



if = ^ 2 {Zf) - eZf = - £(i + mfhal, (7) 

U/i j>0 



(8) 

The relation ([]) demonstrates explicitly the main computational advantage of finding at 
first the polarization function II s (Q 2 ) against a direct calculation of R s (s) in the case 
of massless pQCD. Indeed, in order the derivative gf^n s and, consequently, R s (s) 

depends on the very function IT which is multiplied by at least one factor of a s . Thus, 
one needs only to know IT to order a 2 and the anomalous dimension 7^ (a s ) to order 
to find all Q-dependent terms in II s at 0(af), since the beta function and the quark mass 
anomalous dimension 7m are reliably known to order from Refs. [17-20]. 

As we shall see later both problems are eventually reduceable to calculation of a specific 
class of diagrams representing some massless Feynman integrals (FI) depending on only 
one external momentum (to be named p-integrals below) and with the number of loops 
not exceeding three. Note that this is in an obvious disagreement with a statement from 
Ref. [JD]] about the inevitable necessity to compute finite parts of four-loop p-integrals in 
order to find the 0(a? s ) contribution to R s (s). Yet, we agree with the author of Ref. [11 



that if such a point were correct it would certainly preclude, at least at the present state 
of the art, any possibility of analytical calculation of R s (s) to order a^. 
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3 Calculation of n s 



In order a 2 the polarization operator II s is contributed by twelve three-loop p-integrals. 
Such a calculation is now to be considered as almost trivial one because of three facts: 

a There is an elaborated algorithm which provides a way to evaluate analytically divergent 
as well as finite parts of any three- loop dimensionally regulated p-integral plj |22| . 

b The algorithm is reliably implemented in the language of FORM 
named MINCER in Ref. @. 



as the package 



c The package has been extensively tested and a chance of a bug in it seems to be very 
small. 



We have computed the IIq to a 2 using the general gauge with the gluon propagator 
the following result 



C,^^) I q 2 . On performing the renormalization according to eq. (|3]) we have obtained 



U SS (Q' 



d[R]Q* 



2 hi + a s C F 



131 17 u 2 3, 2 /i 2 



Cl 



i^ +2 4C(3)--C(4) 
64 sv ; 4 SV ; 



... 691 fi 2 
15 C(5) - — In— 



9 u 2 
+ 2^ 



105 2[ x 2 
T6- ln Q^ 



4 Q 2 



32 

2 



+ Cp C £ 



14419 75 A/oN 



| f(4) + | C(5) - 8 '« 



32 ln "Q 2 



31 ., oNl /i 2 71, 2 /i 2 11, 3 M 



+ C F Tn f 



511 



65 



+ 



11 



36- 4C(3) + Y ln ^- 2C(3)ln -Q 2 ' 6-g 2 6-g 2 



111 



2^ 



+ - ln d r 



(9) 



Here a s = a s (fi)/7r, Ca and Cp are the Casimir operators of the adjoint and quark 
(defining) representations of the colour group; T is the normalization of the trace of 
generators of quark representation Tr(t a t b ) = T5 ab ; rif is the number of quark flavours; 
d[R] is the dimension of the quark representation of the colour group. Note, that all the 
Q-dependent terms in (|9|) are in agreement with the results of M. The independence of 
(P) from the gauge parameter £ is of course expected on general grounds. 



4 Calculation of the anomalous dimension 

There is about a hundred of diagrams contributing to II s in order a^. According to (|3]), 
to compute 7^ one should determine the UV counterterms (that is essentially divergent 
parts) of all the corresponding four-loop p-integrals. Unfortunately, at present there exists 
no way to compute directly the divergent part of a generic four-loop p-integral. 



3 



The only available (unfortunately rather involved) approach to perform such calcula- 
tions analytically is to use the method of Infrared Rearrangement (IRR) discovered by 
A. Vladimirov in Ref. |25| (see also Refs. [0]). The method simplifies calculation of 
UV counterterms by the effective use of the following important theorem: in a class of 
minimal renormalization schemes (including MS- and MS-ones) any UV counterterm has 
to be polynomial in external momenta and masses |27J. It amounts to an appropriate 
transformation of the IR structure of FI's by expanding the latter in a formal Taylor 
series with respect to some external momenta and masses, with resulting FI's being much 
simpler to calculate. 

The method of IRR was significantly extended with elaborating a so-called i?*-operation 



in Refs. [28, 22, 30]. By an explicit construction of the corresponding algorithm, it has 



been shown in Ref. p9f how to reduce calculation of the UV counterterm for an arbitrary 



(h+l)-loop FI to evaluation of divergent and finite parts of some appropriately chosen 
h-loop p-integrals. 

It should be stressed that the .Reoperation is absolutely essential for the algorithm to 
work in general case, though in most (but not in all) practical cases one could proceed 
without it. However, such a practice forces the use of a diagram- wise renormalization 
procedure; the latter, being very difficult to perform by a computer, implies a huge amount 
of highly error-prone and time-consuming manipulations with hundreds of diagrams. 

For instance, the only QCD four-loop problem performed by now is the evaluation of 
the ratio R(s) = o"tot( e+e ~ — > hadrons)/<r(e + e~ — > /i + yu~) to order af. It was done in 
Refs. pi| , |3"i2"[ within the diagram-wise renormalization approach. The core of the problem 
is the calculation of the four-loop contribution to the photon anomalous dimension enter- 
ing into the RG equation for the photon polarization operator. The initial 98 four-loop 
diagrams contributing to the photon polarization operator proliferate to about 250 after 
the IRR procedure is applied. In addition these diagrams contain about 600 various sub- 
diagrams which should be computed separately in order to subtract UV sub divergences. 

Technically, the calculation of 7^ s is obviously very similar to that of the photon 
anomalous dimension. Consequently, any attempt of straightforward repetition of those 
calculations for the case of the scalar correlator would mean a few man-years of routine 
and boring work and, thus, would be not acceptable for the present author. 

Below we will use, instead, the power of the -Reoperation and simplify the application 
of IRR to the problem so far that both UV and IR renormalizations can be done in a 
global form and, consequently, can be simply performed by computer. 

We begin with from the Dayson-Schwinger equation for the correlator ([I]) written in 
the bare form[] 

1%,^) = - / dp^Tr[G (p + q ,a s )T s (p, q ,a° s )G (p,a% (10) 

Here G° and are the full quark propagator and the scalar current vertex function 
respectively; below the integration with respect to the loop momentum p with the weight 
function will be only understood but not explicitly displayed. 

1 For simplicity we set the 't Hooft-Veltman unit of mass /x equal to 1 below. 
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The renormalized version of ( |i~0[) reads 



n s (g, a.) = ZfQ 2 - $Tr[G°(p + q, a° s )T° s (p, q, a° s )G°(p, a% (11) 



Here Z2 is the quark wave function renormalization constant; Z$ is the renormalization 
constant of the scalar quark current defined as 

^] = Z s /Z 2 

where the current inside squared brackets is the renormalized one. A well-known equality 
Z$ = Z m Z2 implies the equivalence of (||) and (|TT1). 

From the finiteness of the renormalized correlator one gets 

Zf = -^|^(||n 9 Tr[G (p + g,a°)r 5 (p,g,a°)G (p,a°)] 

+ ^□ gTr[G o (p + g;a o )r o (p)g)a o )G o fea o )] ^| (12) 

where K e f(e) stands for the singular part of the Laurent expansion of /(e) in e near 
e = and 5Z$ — Z$ — 1. In Eq. ([12]) we have let a Dalambertian with respect to the 
external momentum q act on quadratically divergent diagrams to transform them to the 
logarithmically divergent ones. We also have introduced an auxiliary mass dependence 
to a quark propagator — the one entering into the "left" current Js — by making the 
following replacement 

G (p,a° s )^G (p,a s )p 2 /( P 2 -ml). (13) 

Note, please, that the auxiliary mass dependence has caused somewhat more complicated 
structure of UV renormalizations in the right hand side of Eq. ( |T2| ) . 

The idea of the method of IRR is quite simple: since the renormalization constant Zf 
does not depend on anything dimensionful one could significantly simplify its calculation 
by nullifying the momentum q in Eq. ([12]). The only requirement which must be respected 
is the absence of any IR singularities in the resulting integrals. Unfortunately, a mass 
introduced to a propagator distinguished by some topological property like the one we 
created above is not always sufficient to suppress all IR divergences in all diagrams. For 
instance, if q — then there appear completely massless tadpoles in the second term on 
the rhs of Eq. (|12"D . Thus, the eq. ( |i"2|) with q = is not valid unless the unwanted IR 
poles are all identified and subtracted away. 

This is certainly the job the .Reoperation was created for! The rules of Ref. |H| spell 
how to do it on the diagram-wise level. The only remaining problem is to disentangle the 
relevant combinatorics and put down the IR subtractions in a global form. The task is 
facilitated by the fact that, as shown in Ref the IR counterterms for an arbitrary 
diagram can be determined in terms of some properly chosen combination of the UV ones. 
The final formula incorporating all necessary UV and IR subtractions reads 

Zf = -^{^^□ 9 Tr[G (p + g,a°)r 5 (p,g,a°)G°(p,a°)]| g = 




7SS f7 7 7SS ^ 

^ 4 :Tr[^(0 I 0,a^-^^|. (14) 
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Here, by 5T° s (p, q, a°) we denote the vertex function of the scalar current with the tree 
contribution removed. The "tilde" atop S again means that in every diagram the quark 
propagator entering to the vertex J s is softened at small momenta by means of the 
auxiliary mass mo according to Eq. flT3|). The bare coupling constant a° is to be understood 
as ct s = Z a ct s ^ with Z a being the coupling constant renormalization constant. To our 
accuracy in a s 

Z a = l + a s (-/3 /e) + a 2 (/3 2 /e 2 - /3 1 /(2e)) . 

Finally, an inspection of (|T4|) immediately shows that, in order to find the (n + l)-loop 
correction to Z^ s , one needs only to know the renormalization constants Z^ s , Z2 and Z m 
to order well as the bare Green functions 



d_ 



G°( P ,a° s ), — [r° s (p,g,a°)]|_ n , D q [V° s (p, q, a° s )]\„ _ n , 5r|(0,0,a°) (15) 



up to (and including) n-loops. Thus, we have obtained a general formula for Z^ s in terms 
of bare p-integrals with explicitly resolved UV and IR subtractions. 



5 Results and discussion 



We have computed with the program MINCER [24| the unrenormalized three-loop Green 
functions (pl|) as well as the renormalization constants Z m and Z2 to order a^. The 
calculations have been performed in the general covariant gauge. The total calculational 
time with an IBM workstation is about 40 hours for the general gauge; for the Feynman 
one it is reduced to about 4 hours. 

Then we have used Eqs. ([[4]) and (|7|) to find to order a^. Our results for 7^ and 
R s (s) read 



7 SS 



d[R] { -2 



+ CpC, 



+ 



+ 



+ 



a s Cp 
4651 



Cl 



119 9 

32-"2C(3)) + C',C', 



7 ^ + l«Z))+C F Tn f 



384 
/641 _ 259 



f C(3) + f C(4) + f C(5) 



V 48 



16 



C(3) + ^C(4) + f C(5)) 



CpC\ 



267889 475 ., , 33 ^,. 45 ., r . 



+ C|T n/ (^-ic(3)-3C(4) 



/ 6^1 ^ Q 
Cp C A Tn f [ — + - C(3) + - C(4) 
/ V7776 3 SV ; 4 SV ; 

<^(!£-!«»)]}- 



(16) 
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and 



R s (s, n) = d[R] {l + a a (ji) [ Sl + 2^ m ln£ 

+ [s 2 + ln£ (s x /3 + 2 Sl7 ^ + 2 7 ^) + ln 2 ^ (/3 o7 ^ + 2( 7 £J 2 ) 
+ af(//) [s 3 + ln£ (2s 2 /? + «iA + 2s 27 ^ + 2s l7 ^ + 2 7 £J 

+ ln 2 ^ ( Sl /? 2 + 3 Sl f3 0l ° m + (3 l7 ° m + 2 Sl ( 7 ^) 2 + 2MI + HcL) 

+ ln 3 ^(l« + 2/? (7^) 2 + !(7^) 3 )]}. 

(17) 

Here the coefficients of the beta-function /3q, /3i and the quark mass anomalous dimension 

7m) 



7° 



0,1,2 are defined according to (0,^) and read 
1 



[3Ci 



7 1 

/ m 



1 

16 



|c|. + ^C F C A - fc F Tn 



1 

64 



129 



129 



ClC 



H413 _ _ 2 



108 



+C 2 T % (48C(3) - 46) + C F C A Tn f ( -48((3) - — J 



556\ 140 



C F T 2 n 2 



27 



y Ca ~ s Tnf . 



*=16 



34 



20 



— Cl-4C F Tn/-— CUTn, 



At last, the coefficients si, s 2 and S3 are found to be 



(19) 



S3 



17 



+ c% 



+ C F Tn f 

23443 27 



s 2 = C F 
65 1 



-7T 



16 + 12 



691 
.~64~ 

vr 2 + C(3) 



9 



893 11 2 31 



-^ 2 -^C(3) + ^C(5) 
768 16 16 w 8 w 



+ C 2 F Ca 



13153 383 



1089 



-7T 



+ C F C; 
+ C 2 F Tn f 



192 96 
3894493 1715 



32 



C(3) 



145 
Iff 



C(5) 



62208 



864 



-7T 



2329 A/oN 25 , . 



88 65 



65 



46 



-7T 



+ -C(3) + -C(4)-5C(5) 



4 



4 



+ C F C A Tn f 



33475 571 2 22 , . 

+ 7T 2 + — C 3 

432 3 SV ; 



+ C F T 2 n 2 f 



972 
15511 11 



I C(4) + § C(5) 



3888 54 



-7T 



C(3) 



(20) 



We observe that neither 7 ^ s nor R s (s) depend on the gauge fixing parameter £ as it must 
be. For the standard QCD colour group values d[R] = 3, C F = 4/3, Ca = 3 and T = 1/2 
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we get for R s (s) with fi 2 set to s: 



R a {s) = 3<j 1+a 
+a 3 s 



3 



+a 2 s 



10801 19 2 39 65 1 2 2 



17 

" - — vr 2 - — C(3) - — nt + — TT 2 n f + -C(3)n f 

144 12 2 sv ; 24 7 18 i 3 SV ; / 



6163613 3535 2 109735 815 46147 277 2 



-7T 



C(3) + ^ r C(5)--^ r n / + — tHtv 



5184 72 216 " v ' 12 " v 7 486 J 72 

262 5_, 25 ., r , 15511 2 11 2 2 1 ., , 2 

+ C(3) n/ " « C(4) n, - - C(5) n, + — -— vr 2 - - C(3) n) 



(21) 



or, in the numerical form, 



a. As) , , / aufs) 



2 



i? s ( s ) = 3 { 1 + 5.66667^^ + (35.93996 - 1.35865n/) (22) 



7T 




+ (164.13921 - 25.77119/2/ + 0.258974n 2 
At last, for the phenomenologically relevant case of nf = 5 we obtain 

R s ( s ) = 3 J 1 + 5.66667^^ + 29.1467 I J +41.7576(^^1 )> . (23) 



7T 1 \ 7T 



Due to eq. (0) the combination m 2 i? s (s = Mj^) is directly related to the Higgs decay rate 
to the 66 pairs plus gluons. The corresponding expression including power suppressed 
corrections reads 



T(H -> 66) 

3G 

F M H m 2 b (M H )[ 1 + 5.67 a s (M#) + 29.15 a 2 (M H ) +41.76 a 3 s (M H ) 



+ m l^ } (-6 - ^a s (M H ) - 87.72 a 2 (M g 



• (24) 



Let us take as input parameters a A.q OD = 233 MeV and a bottom pole mass of m^ ole = 
4.7 GeV. The latter translates into the running mass m b (M%) = 2.84/2.75/2.69 GeV 
for Higgs masses of Mh = 70/100/130 GeV. All other light quarks are assumed to be 
massless. One arrives at the following values for the strong coupling constant: a s (Mfj) = 
0.125/0.118/0.114 corresponding to the three different values of Mjj. The corresponding 
numerical values for different contributions to R (Mfj) are givenQ in Table 1. 
To summarize: we have suggested a new convenient way to compute the UV renormaliza- 
tion constant of the correlator of the scalar quark currents. Our final formula flli]) directly 



2 It should be noted that this discussion is somewhat inconsistent because of the following reasons. 
First, the relation between the running and the pole masses is known only with the a 2 s accuracy fffi] . 
Second, one needs not yet available coefficients, 7^ and fa, for taking into account the terms of order 
(a s (m^ ole ) — a s (.M#)) 3 in the running of the quark mass and the coupling constant. Nevertheless, the 
numbers in Table 1 are correct for the given values of a s {Mjj) and m^Mfj). 
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a s (M H ) 


m b (M 2 H ) 


(a s (M H )) 


0{al{M H )) 


0(al(M H )) 


70 GeV 


0.125 


2.84 GeV 


0.226 - 0.00262 


0.0461 - 0.00023 


0.0026 


100 GeV 


0.118 


2.75 GeV 


0.213 - 0.00114 


0.0411 - 0.00009 


0.0022 


130 GeV 


0.114 


2.69 GeV 


0.206 - 0.00062 


0.0384 - 0.00005 


0.0020 



Table 1: The contributions of different orders to R s /3; the first number stands for the 
massless part and the second one for the power suppressed 0{^t) part (where available). 

expresses the constant in terms of unrenormalized p-integrals, with all UV and IR sub- 
tractions being implemented in a global form. The formula is ideally suited for carrying 
out completely automatic calculations and can be easily extended for the case of general 
bilinear quark currents. Detailed derivation of the formula will be presented elsewhere. 
Using the formula and the FORM version of the MINCER we have computed the 0(a^) 
correction to the anomalous dimension to the absorptive part of the scalar correlator 
R s (s) and to T(H — > hadrons) in pQCD. All the calculations have been performed with 
the use of the general covariant gauge. The gauge independence of the results constitutes 
a strong check of the correctness of our approach. 

Numerically, the correction of order to T(H — > hadrons) proves to be relatively 
small. However, for the Higgs boson with intermediate mass it is more important than 
the power suppressed 0(-^-a 2 s ) contribution, (see Table 1). 
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